Rosen-Zener Transition in a Nonlinear Two-Level System 
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We study Rosen-Zener transition (RZT) in a nonlinear two-level system in which the level energies 
depend on the occupation of the levels, representing a mean-field type of interaction between the 
particles. We find that the nonlinearity could affect the quantum transition dramatically. At 
certain nonlinearity the 100% population transfer between two levels is observed and found to be 
robust over a very wide range of external parameters. On the other hand, the quantum transition 
could be completely blocked by a strong nonlinearity. In the sudden and adiabatic limits we have 
derived analytical expressions for the transition probability. Numerical explorations are made for a 
wide range of parameters of the general case. Possible applications of our theory to Bose-Einstern 
Condensates (BECs) are discussed. 
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I. INTRODUCTION 

The Rosen-Zener model was firstly proposed to study 
the spin- flip of two- level (Hyperfine Zeeman energy level) 
atoms interacting with a rotating magnetic field by 
N.Rosen and C.Zener to account for the double Stern- 
Gerlach experiments Different from the well-known 
Landau-Zener model that depicts the tunnelling dynam- 
ics between two avoided-crossing energy levels |2|, in the 
Rosen-Zener model the energy bias between two levels 
is fixed and the coupling between two modes is time- 
dependent described by a rectangular Q, Gaussian 
exponential Q or a hyperbolic-secant functions [H . This 
model constantly attracts much attention not only be- 
cause it has exact analytic solution providing a bridge to 
understand complicated multi-mode systems [|| but also 
due to its versatile applications in nonresonant charge ex- 
change of ion-atom collision 0], laser- induced excitation 
[H, nuclear magnetic resonance technique and quantum 
computation[9(, to name only a few. 

In the present paper, we extend the Rosen-Zener model 
to nonlinear case and want to see how nonlinearity af- 
fects the quantum transition dynamics in this system. 
Our work is motivated by recent flourish of Bose Einstein 
condensates (BECs), where nonlinearity naturally arises 
from a mean-field treatment of the interaction between 
particles. In fact, nonlinear effects consta ntly e merge in 
BEC system recently, such as self-trapping [10. 1 1 ll. 1 1 2 . 1 1 3j . 
superfluiditvfliL instabilityflEj and nonlinear Landau- 
Zener tunneling (lH \vf\ . Extending the famous Rosen- 
Zener model to nonlinear case is of great interest and 
worthwhile at this point. 

Our paper is organized as follows. In Sec. II, we intro- 
duce the nonlinear Rosen-Zener model. In Sec. Ill, for the 
degenerate case that the energy bias between two levels 
is zero, our numerical calculations reveal the significant 
effects of nonlinearity on RZT. We then derive analytic 



expressions for RZT in sudden and adiabatic limits. In 
Sec. IV, our discussions are extended to nondegenerate 
case. Numerical explorations for a wide range of pa- 
rameters are made. Interesting phenomena are presented 
and discussed. Finally, in Sec.V, possible experimental 
realization of our model in Bose-Einstern Condensates 
(BECs) are discussed. 



II. NONLINEAR ROSEN-ZENER MODEL 

The nonlinear two-mode system we consider is de- 
scribed by following dimensionless Schrodinger equation, 
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with the Hamiltonian given by 
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where a x and a z are Pauli matrices, 7 and v are the 
energy bias and coupling strength between two modes, 
respectively, c is the nonlinear parameter describing the 
inter-atomic interaction. The total probability |a| 2 + \ b\ 2 
is conserved and set to be 1 without losing generality. 

Different from the Landau-Zener type model where 
coupling keep s constant and energy bias varies in time 
linearly 0, H, [13] , m this model, the energy bias 7 and 
nonlinearity c arc set to be constant whereas the coupling 
v is time-dependent governed by an external pulse field 
of the form [181 ] . 



0, t < 0,t > T 
u sin 2 (f ), t G [0,T] 



(3) 
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where T is the scanning period of the external field. 

In our following study, we assume the quantum state 
is prepared on one mode initially. With the external field 
turned on, quantum transition between different modes 
emerges. What we concern is the population dynam- 
ics in the presence of the external field. The transition 
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probability T is defined as the probability of the particle 
occupying the other mode after the external field pulse 
is over. 

We start our analysis with the simplest case, i.e., 
both the energy bias and nonlinear parameter vanish 
(7 = and c = 0). In this case, the Schrodinger equa- 
tion denoted by Eq.fl]) is solvable analytically. Setting 
(a, b) = (1,0) as the initial condition, we readily obtain 
the probability of the particle populated on the other 
mode as the function of time, 

, , , , s,2 , 2?rf - Tsin^ , , 

p{t)=\b{t)\ 2 =sm 2 {v Q 21), (4) 

then total transition probability is obtained by substitut- 
ing t = T into the above equation, 

r = p(T) = sin 2 (^). (5) 

The above expression demonstrates a perfect Rabi- 
like oscillation of the transition probability versus the 
pulse duration or scanning period T of the external field. 
The oscillation frequency is proportional to the maximum 
coupling strength v$. 

III. NONLINEAR RZT FOR DEGENERATE 

CASE (7 = 0) 

A. General Properties 

With the emergence of nonlinearity, the transition 
dynamics dramatically changes. In this case, the 
Schodinger Eq.JT]) is no longer analytically solvable [lj|. 
We therefore exploit 4-5th Runge-Kutta algorithm to 
trace the quantum evolution numerically and calculate 
the transition probability meanwhile. In our calculation, 
we choose the maximum coupling strength v$ as the en- 
ergy scale, thus the weak nonlinearity and strong nonlin- 
earity refer to c/vq « 1 and c/vq >> 1, respectively. 

Our numerical results are presented in Fig[T] where 
the transition probability as the function of the scan- 
ning period are plotted against the nonlinear parameters 
that range from weak nonlinearity to strong nonlinearity. 
In FigQJa) we see a regular oscillating pattern, agreeing 
with the analytic prediction of formula |5]). This regu- 
lar periodic pattern is destroyed with the emergence of 
nonlinearity. For weak nonlinear case c/vq < 1, the peri- 
odicity is lost only in short pulse regime, i.e., < T < 50 
at c/vq = 0.9. Whereas in the regime of large scanning 
period, a rectangular periodic pattern revives instead. 
The period of the rectangular oscillation increases with 
the enhancement of nonlinearity. The above rectangular 
pattern is of particular interest in practice because it rep- 
resents a 100% population transfer between two modes 
robustly in a wide parameter regime. 

For the case of strong nonlinearity, i.e., c/vq > 1 (see 
Fig. Id), the quantum transition between two modes is 
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FIG. 1: (Color online.) Numerical results for the transition 
probability versus the scanning period under different nonlin- 
ear parameters c/vo (b) 0.6, (c) 0.9, (d) 1.3. For comparison, 
we also include the result of linear case in (a), which well 
reproduces the results predicted by formula J5). 



even more affected by the nonlinearity: The oscillation 
pattern is completely broken, and when the scanning pe- 
riod T > 15 the quantum transition is totally blocked. 

The above phenomena are interesting and intriguing, 
need detailed investigation. Our further analysis includes 
two limit cases, namely, adiabatic limit and sudden limit. 
The adiabatic limit means that the external field varies 
slowly compared with the intrinsic motion of the system. 
From the formula© we see the period of intrinsic mo- 
tion is characterized by Att/vq while the external field 
is characterized by T. Thus, the adiabatic limit means 
T » Att/vq or vqT » Air whereas the sudden limit 
corresponds to vqT << 4ir. In the following sections, 
we will further derive some analytical formulas for the 
transition probability and explain the above phenomena. 



B. Analytic Results for Sudden Limit 

Case(«oT << 47r) 

In the sudden limit that the scanning period of the 
external field is small enough, it is clearly seen that non- 
linearity always suppress the transition from one mode 
to another mode (see Figj2j) . This is because of the com- 
petition between the on-site energy characterized by the 
nonlinear term and the hopping energy characterized by 
the coupling strength. For a certain coupling strength 
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FIG. 2: (Color online). Transition probability in the sudden 
limit for varied nonlinear parameters. The scatters are ob- 
tained by directly integrating Eq.fl}. The solid lines demon- 
strate the analytical results based on Eq.©. They are in 
good agreement. 

v, the larger is the nonlinear interaction c, the larger is 
the on-site energy, therefore it will block the population 
transfer between two modes more strongly and the oscil- 
lation between two modes becomes slower. 

Now we derive some analytic results supporting the 
above argument. Because the transition probability is 
small, an explicit analytic expression can be obtained 
using perturbation theory. We begin with the variable 
transformation, 




As a result, the diagonal terms in the Hamiltonian are 
transformed away, and we have 

i.' n\ f u ° • 2 nt -*/„' c(|&'| 2 -|a'| 2 )di ',, fo , 
b(t) = ysin'-e I ' 1 1 adt. (8) 

The first-order amplitude of b (t) can be obtained by as- 
suming a =1,6 = on the right-hand side of the above 
equation. This yields b' (t) = f* ^ sin 2 ^e ict dt, and the 
transition probability equals to 

27rS 2 [l-cos(cT)] 

c 2 (47r 2 _ C 2 T 2)2 • W 

For linear case of c = 0, the above formula reduces to 
|6(T)| 2 = (^-) 2 , which is the first order expansion of 
the exact solution of transition probability of Eq. ([5|) . For 
the nonlinear cases, this analytical approximation is com- 
pared with our numerical results in FigJ5]and good agree- 
ment is shown. Expanding Eq.® using c as the small 
parameter gives, 

r^d-fff!,, (10) 

which indicates that the quantitative decrease of transi- 
tion probability is proportional to the square of nonlinear 
parameter c. 
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FIG. 3: Typical eigenenergy level structure in (a) linear case, 
(b) weak nonlinear case, (c) strong nonlinear case. The 
dashed lines correspond to unstable eigenstates [2p| ]. 



C. Analytic Results for Adiabatic Limit 

Case(i;oT >> 4n) 

According to the adiabatic theory of nonlinear quan- 
tum mechanics [21] , the characters of quantum transition 
in the adiabatic limit should be entirely determined by 
the structure of the energy levels and the properties of 
corresponding eigenstates. The eigenstates of the system 
satisfy that, 




Solving the above nonlinear equations together with total 
particle conservation condition \a\ +\b\ =1, we readily 
obtain the chemical potential [i and the eigenstate (a, b). 
The eigenenergies can be derived according to the rela- 
tionship e = fi — c/2(\a\ + |6| 4 )- Their dependence on 
the parameters is plotted in Fig[3] for the cases of lin- 
earity, weak nonlinearity and strong nonlinearity, respec- 
tively. Striking phenomena are induced by the nonlinear- 
ity: Firstly, the structure of the energy levels change dra- 
matically. In the linear case, there are two energy levels 
that are symmetric about a horizontal axis (see Fig. 3a). 
However, the symmetry breaks down in the presence of 
nonlinearity and a new branch of eigenenergies emerges. 
For the weak nonlinearity (i.e., Fig. 3b), at two ends and 
near the peak of the field pulse there exists two levels, in 
other regime there are three energy levels. When the non- 
linearity is strong (i.e., Fig. 3c), apart from the two ends, 
the number of the energy levels are three. Secondly, the 
eigenstates of the mid level ( e.g. denoted by the dashed 
line in Fig[3]) are unstable. This is evaluated by investi- 
gating the Hamiltonian- Jaccobi matrix obtained by lin- 
earizing the nonlinear equation (|11|) around the eigen- 
states. The eigenvalues of the Hamiltonian- Jaccobi ma- 
trix can be real, complex or pure imaginary. Only pure 
imaginary eigenvalues correspond to stable states, others 
indicate the unstable ones [2Cj. 

The above changes of the energy level in the presence of 
nonlinearity is expected to affect the quantum transition 
between levels greatly. However, because of the degen- 
eracy of our concerned states of (1,0) and (0,1), from 
the above diagram of the energy levels we cannot draw a 
definite conclusion about the adiabatic evolution of the 
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dH/ds,ds/dt= -dH/86, 
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FIG. 4: Evolution of the phase space motions as c/v changes 
adiabatically (a)10, (b)2, (c)l, (d)0.5. The second column is 
the corresponding energy curve at 6 — (thin line) and tt 
(heavy line) . The arrows indicate the shifting direction of the 
fixed points as v increases. 



initial state (1,0). In the following, we study an equiv- 
alent classical Hamiltonian instead, and with analyzing 
its phase space we achieve insight into the adiabatic evo- 
lution of the above nonlinear quantum system. 

Following the theoretical methodology proposed in 
Ref. [2l[ we construct the effective classical Hamiltonian 
with introducing two quantities: p — \b(t)\ as the proba- 
bility of particles staying on (0,1) mode , and Q = Q a — Qb 
as the relative phase of the two modes. They form a 
pair of canonical variables of following classical Hamil- 
tonian and satisfy the canonical equations, i.e., dO/dt = 



H 



Vl>(l-p)cosff- 7 (2p-l) s 



(12) 



The above classical system is capable to totally describe 
the dynamic properties of nonlinear quantum Rosen- 
Zener system (1) on a projective Hilbert spacefill . Its 
fixed points, i.e., energy extrema of the classical Hamil- 
tonian, correspond to the quantum eigenstates. For ex- 
ample, in FigHJ the stable elliptic fixed point Pi corre- 
sponds to the upper level of Fig[3J the energy of Pi and 
Pi are identical, therefore they correspond to the same 
energy level, i.e., the lower one in Fig|3j the saddle point 
P3 is unstable, corresponding to the mid level denoted 
by a dashed curve in FigJ3] The adiabatic evolution of 
the quantum eigenstates therefore can be evaluated by 
tracing the shift of the classical fixed points in the phase 
space when the parameter v varies in time slowly [171]. 

The analytic expressions of the fixed points are ob- 
tained from p = and 9 = 0: 



0,tt; 



(13) 



v(l - 2p*) 



cos0*+c(l-2p*)=O. (14) 



The number of fixed points depends on the instantaneous 
coupling strength v and the nonlinear parameter c. If 
c/v < 1, there exists only two fixed points (9*,p*) = 
(0, 1/2), (?r, 1/2). However, when c/v > 1, there exists 
four fixed points (Pi, P 2 , P3 and P 4 in Figf4|a)(b) ) : 

(6*,p*) = (0, 1/2), (tt, 1/2), (tt, 1(1 ± y/l-£)). One of 

them (P3) is a saddle point, while the other three (Pi, P2 
and P4) are all elliptic points corresponding to the local 
maximum (Pi) and minimum (P2 and P4) of the classical 
Hamiltonian. 

When we increase v, P2, P3 and P4 will merge into a 
new stable fixed point P234 in a regime satisfying the con- 
dition v/c > 1, as can be seen from FigQJa^b-^c-^d) . 
An interesting question arises when the parameter v de- 
creases to below c again: Which point will the state 
choose to follow when P234 bifurcates into P2, P3 and 
P 4 (FiggH^c^b^a). The state that totally follows P 2 
will give a zero of the adiabatic transition probability, 
whereas the state that totally follows the P4 will cor- 
respond a complete population transfer. This classical 
picture could explain why we see a rectangular pattern 
in FigEEJ 

The above scenario is further supported by directly 
tracing the evolution of particles as shown in FigEJ where 
we also demonstrate the temporal evolution of the fixed 
points. In the early stage, we see that the state firmly 
follows the fixed point P2. It starts to show a small oscil- 
lation when the fixed points P2, P3 and P4 merge. After 
that, the state either follows the fixed point P2 or P4. The 
interesting thing is that, at certain parameter, a slight 
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2000 



FIG. 5: Evolution of the particles and the fixed points at 
c/vo = 0.9. Entirely different transition probabilities are 
observed for a slight variation of the scanning period: (a) 
T=2024, (b) T=2025. 



change on the period T could thoroughly change the fi- 
nal transition probability, a signature of the appearance 
of bistable state (Fig©' ->b'). 

Our study shows that whether the state follows P 2 
or P4 is determined by the total dynamical phase accu- 
mulated during the oscillation motion around P234, i-e. 
from t* to t**. Here, t* is the time when P 2 , P3 and P4 
merge into P234, and t** is the moment when P234 bifur- 
cates into P2, P3 and P4. They are obtained by setting 

c = v(t*) = v(t**y. 



T 



t = —sin . — , 

t** = T-t* 



(15) 



To obtain the total phase, we need to calculate the 
instantaneous frequency that characterize the oscillations 
around the fixed point at first. To this end, we expand 
the classical Hamiltonian around the fixed point with p = 
1/2 + S p , 9 = 7r + 5g , leading to 



(16) 



with ignoring the higher order terms. The instantaneous 
frequency is then derived as 



1 



<»(*)= 2V«(*)K*)-c] 

Integrating ui(t) from t* to t** gives the total phase 



(17) 



<p - 



uj{t)dt 



vo - c 



T. 



(18) 
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FIG. 6: Periods of rectangular oscillation as shown in Fig[T] 
under adiabatic limit. Numerical results are obtained by di- 
rectly observing Fig[T]while the analytical curve is the plot of 
function 4n/(vo — c), here vq = 1. They agree well with each 
other. 



This expression indicates that the total phase increases 
linearly with the scanning period. A 7r-value change on 
the phase will change the choice of the state that either 
follows fixed point P 2 or P 4 . Thus the period of rectan- 
gular oscillation observed in Fig[T] under adiabatic limit 
can be expressed as 



^rectangular 47t/(vq c). 



(19) 



To check the above theory, we have numerically solve 
the nonlinear Schrodinger equation for a wide range of 
the parameters. The comparison between the analytical 
result and the numerical data shows a good agreement 
in FiglHl 

The complete suppression of quantum transition un- 
der adiabatic limit in Fig[TJd) can also be explained from 
the above picture. We briefly state it as follows. For the 
strong nonlincarity that c/vq > 1, the phase space evo- 
lution only undergoes FigQJa^b— >a) as v increases and 
decreases. During the process, no collision between the 
fixed points occurs. Thus the state initially populated on 
P2 can safely remain on the fixed point, and finally come 
back to the original state smoothly. Thus no transition 
is observed. 



IV. NONLINEAR RZT FOR 
NONDEGENERATE CASES (7 / 0) 

In this section, we extend our discussions to the nonde- 
generate case (7 / 0). The transition probability is con- 
tour plotted with respect to the energy bias and scanning 
period for different nonlinear parameters ranging from 
weak nonlinearity to strong nonlinearity in Fig[7j 

The plotting reveals main features of the RZT in the 
nondegenerate case: i) From the first panel of Figj7l the 
contour plot shows fringe structure: The bright zones cor- 
respond to high transition probability whereas dark areas 
indicate low transition probability. The fringe pattern is 
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FIG. 7: (color online.) The contour plot of RZT probability 
as the function of energy bias and scanning period at different 
nonlinear parameters (upper panel) and its adiabatic counter- 
part (lower panel). Initially, all the particles are prepared on 
the high energy mode (1,0) (left column) or low energy mode 
(0, 1) (right column), respectively, before the external field is 
turned on. 



formed because of the periodic behavior of the RZT with 
respect to the scanning period. For linear (i.e., Fig[7K) 
and weak nonlinear cases (i.e., Fig[7J;), this periodicity is 
broken by a larger energy bias; For strong nonlinear case 
(i.e., Fig[7k), the periodicity is absent both for larger and 
smaller energy bias. ii)The nonzero energy bias leads to 
the asymmetry of RZT plottings for the initial states of 
(1,0) and (0,1). This is clearly seen by comparing the 




FIG. 8: Energy curve of H ( Eq.fEfl ) at c = 1,7 = 0.1, 
= (thin line) or tt (heavy line), v = (a)0.1, (b)0.5, (c)1.0. 
With the increase of v, P3 and P4 collide with each other and 
disappear subsequently, while P2 maintains during the whole 
process. 



right column with the left column of Fig. 6. iii) The sec- 
ond panel of Fig. 6 shows the contour plot of transition 
probability in the adiabatic limit i.e., very large scan- 
ning period T . In this situation, the fringe structure is 
confined in a zonary region. This zonary region is of 
particular interest in practice because only in this re- 
gion a robust population transfer between two modes is 
possible. Outside the zonary region, the transition prob- 
ability almost equals to zero. Interestingly, the area of 
the zonary region broaden with increasing the nonlinear 
parameter (Fig[7^' — >c' — >e'). When the nonlinearity is 
strong (i.e., Fig.6e'), the zonary region shift away from 
the origin. That is to say, to realize robust population 
transfer for the highly nonlinear case, it requires a large 
energy bias between two modes to compensate the differ- 
ence of the on-site energy caused by the nonlinearity. iv) 
In the zonary region, for a fixed 7, the transition proba- 
bility oscillates rectangularly versus the scanning period 
similar to the behavior shown in Fig.l. But the am- 
plitude of the oscillation decreases monotonously with 
increasing the energy bias 7. This is indicated by the 
gradual darkness from left to right side in the figures. 
The correspondence between oscillating amplitude and 
energy bias suggests that the probability of the popula- 
tion transferred to the other mode can be designed at 
will by tuning the energy bias. 

We can achieve insight into the above findings by an- 
alyzing the evolution of the fixed points corresponding 
to energy extrema. The energy curve in Fig[S] is tilted 
by the energy bias (for comparison one can recall Fig[?]), 
directly leading to the asymmetry on the diagram of tran- 
sition probability. The reason is stated as follows. Based 
on our pictures in Sec.III(C), in the adiabatic limit, the 
transition between two modes is always accompanied by a 
collision between fixed points. However, with increasing 
v, the fixed point corresponding to the low energy mode 
(P2) do not collide with any other fixed points, thus there 
is no transition from low energy mode to higher one in 
the adiabatic limit. The situation is quite different for 
the high energy mode (P4), which collides with the un- 
stable state (P3) and disappear subsequently as shown in 
FigJS^b). The condition for the occurrence of the collision 
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is given by [T3|, 

u=(ci-7§)i (20) 

Considering v can only varies in the interval [0, vq], we 
therefore determine analytically the boundaries of the 

zonary regions as, [0, c] if c < vo and [(c§ — Vq )',c] if 
c > vq. These two analysis is in a good agreement with 
the numerical results ( see Fig|8fc')(e') ). 

Our above discussions assume that the nonlinear pa- 
rameters are positive that corresponds to the repulsive 
interaction between particles. For the attractive inter- 
action, the transition from (1,0) mode to (0,1) mode is 
equivalent to that transition from (0,1) mode to (1,0) 
mode for the repulsive interaction. Therefore, our discus- 
sions can be readily extended to the negative-nonlinearity 
case that corresponds to attractive interaction between 
particles. The detailed discussion is not repeated. 

V. APPLICATIONS AND DISCUSSIONS 

Our model is applicable to describe the quantum tran- 
sition of two-mode BECs. One example is the BECs 
in a double-well. In such system, the wavefunction can 
be described by a superposition of two states that local- 
ize in each well separately. Then, the transition dynam- 
ics of BECs between the two wells is well described by 
Eq.(l)[ll|. The optical double- well can be created, for 
example, by superimposing a blue-detuned laser beam 
upon the center of the magnetic trap [22] . In this case, 7 
denotes the difference of the zero-point energy between 
two wells, c is the interaction between BEC atoms that 
can be adjusted flexibly by Feshbach resonance, v repre- 
sents the height of the barrier that separate the two wells. 



The barrier height can be effectively controlled by adjust- 
ing the intensity of the blue-detuned laser beam. Initially, 
we upload BECs atoms into one well, then ramp up and 
down the barrier slowly, the nonlinear Rosen- Zener tran- 
sition should be observed. 

Another promising candidate to observe nonlinear 
Rosen-Zener transition is a spinor BEC in an optical 
trap. In such case, a near-resonant field is applied to the 
condensates to couple two hyperfine states of 87 Rb, e.g ., 
\F =l,m F = -1) and \F = 2,m F = +1) like in Ref . Hi . 
Within the standard rotating wave approximation, the 
system Hamiltonian can be cast into the similar form 
of Eq.© [2J], where 7 denotes the frequency detun- 
ing, c comes from the contributions of both the homo- 
species and inter-species s-wave scattering, v is propor- 
tional to the intensity of the near-resonant laser field. Af- 
ter preparing the BEC atoms in one internal state, a sin 2 
enveloped laser pulse is shined on the condensates, the 
nonlinear Rosen-Zener transition is expected to emerge. 

In conclusion, we have investigated the Rosen-Zener 
transition in a nonlinear two-level system and show that 
the nonlinearity could dramatically affect the transition 
dynamics leading to many interesting phenomena. The 
nonlinear RZT is suggested to be observed in the two- 
mode BEC systems. 
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